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Abstract - Run-and-tumble dynamics is a wide-spread mechanism of swimming bacteria. The 
accumulation of run-and-tumble microswimmers near impermeable surfaces is studied theoretically 
and numerically in the low-density limit in two and three spatial dimensions. Both uni-modal and 
exponential distributions of the run lengths are considered. Constant run lengths lead to peaks 
and depletions regions in the density distribution of particles near the surface, in contrast to 
exponentially-distributed run lengths. Finally, we present a universal accumulation law for large 
channel widths, which applies not only to run-and-tumble swimmers, but also to many other kinds 
of self-propelled particles. 


Introduction. — Swimming bacteria like E. coli and 
Salmonella, with a body length of just a few micrometers, 
are too small for spatial sensing of a stimulus gradient 
along their body size 00 - Therefore, they have to re¬ 
sort to temporal sensing, where the gradient is determined 
along the swimming trajectory. These bacteria have de¬ 
veloped a procedure of intriguing simplicity for chemotac- 
tic motion — they perform a run-and-tumble motion, in 
which nearly straight swimming segments are interrupted 
by tumbling events, where the run length then depends 
on the sign of the stimulus gradient . There is an in¬ 
teresting connection of this run-and-tumble dynamics to 
Levy flights [^[^, which suggests that this process could 
be a very efficient search strategy [^. 

For a dilute suspension of microswimmers in a bulk 
fluid, run-and-tumble dynamics is strictly equivalent to 
passive-particle diffusion for long times. Here, the diffu¬ 
sion coefficient is given by D^ff = v'^ /tt, where v is the 
swimming velocity, and is the run time [^. The effec¬ 
tive diffusion coefficient is typically much larger than the 
thermal diffusion coefficient. The equivalence also holds in 
the presence of a slowly varying external potential. This 
implies that active Brownian particles (ABPs), which dis¬ 
play a rotational diffusion instead of tumbling events, are 
equivalent to run-and-tumble particles (RTFs) under these 
conditions. 

In fact, the equivalence of ABPs and RTPs has been 


discovered much earlier for the mathematically equivalent 
case of the conformations of semi-flexible polymers. Here, 
the worm-like chain model corresponds to ABPs, whereas 
the freely-jointed chain model corresponds to RTPs. The 
equivalence of the two is expressed by the Kuhn length 
which is the segment length of the freely-jointed chain, 
to equal twice the persistence length of the worm-like 
chain, such that the end-to-end distance is the same in 
both descriptions [^. 


At higher densities of microswimmers, a density- 
dependent motility can cause phase separation and ac¬ 
cumulation of both RTPs [9 -11 and ABPs 12 ■ 14 , which 
indicates that the dynamics of active particles is no longer 
equivalent to passive diffusion. Asymmetric potentials can 
cause rectification of bacterial motion [7 15, fl^. Also, 


walls and obstacles break the diffusion equivalence, be¬ 
cause microswimmers accumulate at walls, in contrast to 
passive particles. Explanations of this surface trapping 
usually invoke hydrodynamics 17 Whether it is the 


detailed hydrodynamics of the corkscrew motion of E. coli 


flagella 19 or the snake-like motion of the sperm tail 20 


or the far-held hydrodynamics of a hydrodynamic dipole 
17 , hydrodynamics provides an effective attraction to¬ 
ward boundaries 18 21 22 . However, for E. coli, noise 
also plays an important role and may even dominate over 
the rather weak hydrodynamic interactions 23 . Further¬ 
more, it has been shown that persistent motion drives 
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swimmers to the wall, even in presence of strong orien¬ 
tational fluctuations [24}p^ . For harmonic confinement, 
accumulation away from the center has also been found for 
run-and-tumble particles [^. Thus, it is not obvious under 
which conditions the equivalence between RTFs, ATPs, 
and passively diffusing particles holds. We want to clarify 
this question from the point of view of wall accumulation 
and confinement. 

In this letter, we investigate the effect of confinement for 
particles with a pure run-and-tumble dynamics in the low- 
density limit and in the absence of hydrodynamic interac¬ 
tions, both analytically and numerically. The structure of 
the density patterns of RTFs at hard walls is found to de¬ 
pend strongly on the dimensionality of the accessible space 
— between two planar walls in three dimensions (3D), or 
along a surface with lateral confinement in two dimensions 
(2D) — and on the run-length distribution, both quanti¬ 
tatively and qualitatively. Here, the relevant parameter is 
the dimensionless ratio between channel width and (aver¬ 
age) run length, whereas propulsion velocity and tumbling 
frequency only enter indirectly via the run length. RTFs 
are predicted to behave quite differently from ABFs. For 
narrow channels and constant run lengths, the distribution 
of tumbling events develops pronounced extrema, with a 
depletion layer near the wall and a maximum at larger 
distances determined by the run length. These structures 
disappear for exponentially distributed run lengths. Thus, 
the behavior depends sensitively on the run-length dis¬ 
tribution. In contrast, for wide channels, we predict a 
(nearly) universal wall-accumulation law for self-propelled 
particles. This wall-accumulation law only relies on sym¬ 
metries and dimensional arguments, and thus holds for 
many different types of microscopic swimmers. 

Wild-type E. coli have an average run length of 12^m 
[30| , which is the same order of magnitude as the chan¬ 
nel width of microfluidic devices used to manipulate and 
study these bacteria 31-33 . Therefore, our results are 
relevant, inter alia, for the design of microfluidic devices 
for rectification and sorting of run-and-tumble bacteria. 


Model and Simulation Technique. — We study 
run-and-tumble dynamics of individual microswimmers in 
confinement. A particle performs a forward run with a 
velocity v for a time r^. Fach run is followed by a tum¬ 
ble event, where a new orientation angle 9 (see Fig. is 
chosen randomly on the unit circle (2D) or unit sphere 
(3D), i.e. there is no memory of the orientation before the 
tumbling event. The particle coordinate z perpendicular 
to the wall is then updated by z{t + Tr) = z{t) -\- cos{6)L. 
Here, the run length L = vTr is either constant, or drawn 
from an exponential distribution depending on the dynam¬ 
ics studied. It is important to note that properties of RTFs 
do not depend on v and separately, but only on the run 
length. Due to symmetry, motion parallel to the wall does 
not have to be considered. If the particle hits a wall, it 
remains there - possibly sliding parallel to the wall - until 
the next tumbling event occurs. After a sufficiently long 



Fig. 1: Schematic of run-and-tumble dynamics. A “run” of 
the particle with a velocity v for a time is followed by a 
tumbling event, resulting in a new orientation 0 . The particle 
is confined between two parallel walls a.t z = ±d. 2 * denotes 
the distance from the walls. 


equilibration time, the probability density to find the par¬ 
ticle at a position z is recorded by a histogram over 10® 
to 10® tumbling steps. A few examples of density distri¬ 
butions for various run lengths are shown in Fig. 

The trajectory of an RTF is completely defined by the 
location of the tumbling events, since the motion between 
these events is just ballistic. This implies, in particular, 
that no orientation vector of the particle is needed to de¬ 
scribe the dynamics. Thus, the continuous-time dynamics 
of a RTF in three spatial and two orientational dimensions 
is mapped onto a discrete-time-step model in one spatial 
dimension. Fhysically and mathematically, the fundamen¬ 
tal quantity to compute in the steady state is then the 
tumbling density 4>{z). The particle density p{z) then fol¬ 
lows from (/>(z) by a convolution, as explained in detail 
below. Both densities are directly asccessible experimen¬ 
tally; however, the tumbling density is more difficult to 
measure, because the trajectories of (all) particles have to 
be traced. 

Thus, we first focus on the more fundamental tumbling 
density, which is the (normalized) probability to find a 
tumbling event at a position z. The mirror symmetry of 
the system is reflected in the symmetry of the tumbling 
density, (j){z) = (l){—z). The time evolution of the tumbling 
density is determined by 

p-\-d 

(j){z,t + Tr)= / (j){z',t)p{z — z')dz'. (1) 

J-d 

Here, p{ Az) is the transfer function of particles moving to 
a new position, which depends implicitly on the run-length 
distribution. It is the number of orientational microstates 
of an unconfined particle which are compatibe with a given 
Az-displacement. This probability density depends on the 
dimensionality of the system, and on the run length dis¬ 
tribution (unimodal or exponentially distributed). At the 
walls, particles accumulate in a J-distribution because all 
particles that hit the wall are located there. Thus, we have 
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Fig. 2: Particle density distribution p{z) for various run 

lengths L. The solid line is the exact analytic solution ( |19[ ) for 
the particle density in three dimensions for run lengths larger 
than the channel width 2d. 


the boundary conditions 

(j){±d,t) =0.5(j)sit)6{z ± d) ( 2 ) 

nd 

0.50s(t-h Tj.) = / ,t)P{—z'— d)dz' (3) 

J-d 


where 0s is the probability to hnd a tumbling event at the 
wall, and P is the cumulative distribution function of p, 
i.e. P{z)=r_^p{z')dz'. 

Constant Run Lengths. — We begin our analysis in 
the simplest case of constant run length L. The transfer 
function is discontinuous at the run length, as runs longer 
than L cannot occur. In three dimensions, the transfer 
function p{z)dz is obtained by an integral over the surface 
of a sphere of radius L with values of the vertical displace¬ 
ment between z and z + dz. This yields immediately that 
the transfer function is 

P(3,c){^) = z)e{L +z), (4) 


where 0(z) is the Heaviside step function. Subscripts de¬ 
note dimensionality and the type of run length distribu¬ 
tion (c=constant, e=exponential). Similarly, in two di¬ 
mensions, integration over a circle of radius L with dis¬ 
placement between z and z -|- dz yields even a divergence 
of the transfer funtion at the run length, 


P(2,c){z) = 


ttLa/I - (z/L)'^ 


Q{L — z)Q{L + z). (5) 


The simplicity of the 3D transfer function allows for an 
analytic solution for narrow channels with 2d < L, 


0(3.c)(^) = ^ + ^-^[5{z -d) + 5{z + d)] (6) 


In 2D, an analytical solution can only be obtained by as¬ 
suming that the number of tumbling events in the bulk 
is negligible compared to tumbling events at the wall (i.e. 


S-d+e4'{z-)dz <C 0s with an infinitesimal length e). This 
results in 


<l>( 2 ,c){z) = 


TT-l-sin ^(2d/L) 

I 


[0.5<5(z — d) -I- 0.5(5(z -I- d)] 


27r-|-27rsin ^(2d/L) 
1 


-b 


_ rf )2 ^ L 2 _ + rf )2 


■ (7) 


These analytical results and corresponding simulation 
data are displayed in Fig. The comparison shows that 
the solution ([^ in 3D and the approximate expression Q 
in 2D work very well for the appropriate regimes. Fig¬ 
ure reveals that the walls induce a very rich structure 
of the tumbling density in the channel for 2d > L, i.e. 
for channels wider than the run length. The density pro¬ 
files all collapse onto a single master curve when the tum¬ 
bling density is scaled by the bulk density 0b (the density 
far away from the wall) and distances are scaled by the 
run length. In this case, the high particle density at the 
walls generates depletion regions near the walls, and two 
pronounced peaks at a distance L from the wall, which 
can easily be recognized in Fig. (middle) and (right). 
In 2D, these primary peaks generate secondary peaks for 
d > L, which are again displaced by a distance L further 
away from wall they first came from. In 3D, the primary 
singularities are too weak to generate visible secondary 
peaks. All bulk singularities disappear for L > 2d, com¬ 
pare Eq. (§, because particles can move directly from one 
wall to the other in a single step. These depletion zones 
and peaks can be understood by starting from a uniform 
bulk distribution plus d-peaks at the walls, and iterating 
Eq. 0 once, as explained in more detail below. 

Exponential Run-Length Distribution. — In the 

case of a distribution of run lengths, the transfer functions 
are obtained by convolution of the (conditional) probabil¬ 
ities p{z\L') — resulting from step with run length L' — 
with the run-length distribution pien{L'), i.e. 

poo 

P{n,len){z^ = / ^(n,c){z\L )p/en(L )dL (8) 

J z 

The integral has a lower boundary at z, because Az = z 
cannot be achieved with L' < z. 

We focus here on exponential run length distributions, 


P/en(7v') = ^exp(-AL'), (9) 


with {L') = L = I/A, which mimic the run-length distri¬ 
bution of E. coli. This yields the transfer functions 


P{ 2 ,e){z) 

P{3,e)[z) =^Ei 



( 10 ) 


P-3 



















Jens Elgeti and Gerhard Gompper 




0.001 0.01 0.1 1 10 100 
d/L 


1.4 

1.2 

3- 

^ 1 
0.8 
0.6 

1.4 

1.2 

3- 

^ 1 
0.8 
0.6 


i 2D 

1 i 

Exponential / 
/ 


d/L=2.5 
d/L=5 
d/L=10 
d/L=2.5 
d/L=5 
d/L=10 . 




' ■ 

' Constant 


0.5 


1 . 1.5 

z/L 


3D 


tu Exponential 

\ “u z ■ I 

\ “s. Z.- I 

\ ■ I 


d/L=2.5 

d/L=5 

d/L=10 

d/L=2.5 

d/L=5 

d/L=10 


0.5 


1 . 1.5 

z/L 


2.5 


2.5 



Fig. 3: Tumbling density profiles (j){z) for (top) two and (bottom) three dimensions, (left) Scaled surface density 4>sd/L as 
a function of the ratio djL oi channel width and run length. Solid lines are analytical approximations for narrow channels (see 
text), dashed lines are fits to the large-channel approximation | |14[ ). Note that for d ^ L, (ipsd/L) approaches a/2, (center) 
Scaled tumbling density (f){z*)/(j)b as a function of the scaled distance z* /L = {z + d)/L from the wall. Dashed lines show the 
approximation which is obtained from a hrst-order iteration of Eq. 0 (see text), (right) Density of tumbling events 

(j){z/L) inside the channel for various (constant) run lengths. For 2d > L, the presence of the two walls induces strong density 
modulations. 


with Kq a Bessel function and Ei an exponential-integral 
function. The continuity and strong decay of these trans¬ 
fer functions leads to the disappearance of all the singular¬ 
ities of the tumbling density found for constant run length, 
see Fig. Thus, the tumbling density is highly sensitive 
to the run-length distribution. 

We estimate the density profile in small channels, by 
neglecting the particles in the bulk in Eq. 0, as for con¬ 
stant run length. In two and three spatial dimensions, we 
then obtain 


Scaling Behavior of Wall Density. — For channels 
much wider than the (average) run length, we can use 
scaling arguments to determine the wall accumulation of 
particles. For L, the tumbling density profile eventu¬ 
ally becomes flat far from the walls. Everything else fixed, 
the surface density has to be linear in the bulk density (j)f, 
defined as the (constant) density far from the walls. Since 
the (average) run length L is the only relevant length scale 
near the wall, the proportionality factor has to be linear 
in L, so that0 


^s(2,e) — [1 XdKQ{2Xd)L—i{2Xd) 


4’s{3,e) 


+XdKi{2Xd)Lo{2Xd)]~^ 
^ l-F2(2Ad) r^ 


( 11 ) 

( 12 ) 


where the Li{x) with i = —1,0 are Struve functions, and 
Ei(x) with 1 = 1,2 are exponential-integral functions. Fig¬ 
ure 1^ demonstrates that these approximations work very 
well for narrow channels. However, the critical length 
scale, where they break down, is significantly lower than 
for constant run length. The reason is that the transfer 
functions for exponentially distributed run lengths decay 
much faster than the transfer functions for constant run 
lengths. This implies that there are more particles inside 
the channel, so that the approximation (f>{z)dz <gC (ps 

breaks down already at smaller channel widths. 


(ps = aLpb 


(13) 


with a dimensionless prefactor a. In a channel of finite 
with d, normalization then gives 


4^s 


aL 

aL + 2d 


1 

1 -b 2d/{aL) ■ 


(14) 


The surface accumulation factor a is independent of run 
length, and only depends on the dynamics (i.e. 2D/3D, 

^Alternatively, it can be argued that the surface density 0s, which 
is dimensionless in our description, can only depend on the ratio of 
the two available length scales L and d, which implies (j)s = F{d/L)^ 
with some unknown scaling function F{x). 

This can also be seen by considering an infinite half-space, with 
a wall at 2 = 0. In this case, the boundary condition is that the 
density approaches 0^ for 2 ^ oo. Then, L is the only available 
length scale. For finite but very wide channel, the density profile 
should not change. However, the normalization of the probability 
density introduces a constraint on 0^, which implies 0^ ~ 1/d. 
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a 

2D 

3D 

RTFs, constant L 

RTFs, exponential L' distr. 

1.01(2) 

1.40(2) 

0.82(2) 

1.14(2) 

ABPs, L = 2^p 

0.80(3) 

0.37(3) 


Table 1: Accumulation factor a for various self-propelled par¬ 
ticles in two and three spatial dimensions. For ABPs, a values 
are obtained from direct Langevin simulations; in this case, we 
employ the “Kuhn” length as the characteristic length scale. 


constant run length/ exponential run-length distribution). 
From our simulations, we obtain the accumulation factors 
shown in Table The large-distance approximation ( |l4| ) 
works very well for d > L, and even for smaller channels 
it is not too far off (see Fig. [^. Unimodal and exponen¬ 
tial run-length distributions result in accumulation factors 
a, which are clearly different, but still of the same order 
of magnitude. Thus, measurements of a = (ps/iL(j){,) for 
d ^ L might provide a new possibility to characterize 
run-length distributions experimentally. 

Since these arguments rely only on dimensional anal¬ 
ysis, the results should be valid for other types of self- 
propelled particles as well (as long as there is one domi¬ 
nant length scale of the dynamics). For active Brownian 
spheres, this length scale is the persistence length of the 
trajectory = vjDr, where Dr is the rotational diffusion 
coefficient. Using data from Ref. [^, we find indeed an 
excellent agreement for channels much larger than the dif¬ 
fusive length scale Id = sjDtjDr, where Dt is the transla¬ 
tional diffusion coefficient (see Fig. |^. Note that Eq. (14) 
also predicts a crossover from narrow- to wide-channel be¬ 
havior aX L ~ 2d for all kinds of self-propelled particles. 
The fact that the a-values in Table differ significantly 
for RTFs and ABPs clearly demonstrates that these two 
types of self-propelled motion are not equivalent near sur¬ 
faces. However, the fact that these factors are all of order 
unity emphasizes the generic aspect of wall accumulation. 

Near-Wall Density in Wide Channels. — For 

wider channels, the surface density c/g is well described 
by Eq. To understand the structure of the density dis¬ 
tribution close to the wall, the stationary form of Eqs. (1) 
to (3) can be used to obtain an analytical approximation. 
We start as an initial guess with a ^-distribution at the 
wall, with an amplitude 4>s, plus a constant tumbling den¬ 
sity (ftb in the bulk (see Eq. 14). An iteration with Eq. (1) 
then yields 

= ^6{z*) + 4>b6{z*) [ P{z')dz' 

4 Jo 

+ yP(^*) j ~ ■ (15) 



d/^p 


Fig. 4: Surface excess s = f (p(z) — pb)dz for self-propelled 


Brownian spheres. Data from Ref. 
sidered equivalent to 4>s- The wide-channel approximation (141 


Here, s is to be con- 
el approximation (141 
works already rather well for channels widths 2d comparable 
to the persistence length ^p, as long as ^p is much larger than 
the diffusive length scale Id = \/DtjDr- 


ical results for the structure of the tumbling-density pro¬ 
file. In particular, for constant run lengths, it reproduces 
and explains the near-wall dip in the tumbling density, 
i.e. the formation of a depletion layer close to the wall, 
and in turn a peak and discontinuity at z = L. As shown 
in Fig. 1^ (right), this peak leads to interesting patterns 
in the tumbling-density distribution for channel widths 
larger than the run length, in particular a very pronounced 
peak in the channel center for d = L. 

Microswimmer Density. — Finally, we connect the 
tumbling density to the number density of microswim¬ 
mers. This requires the convolution of the tumbling den¬ 
sity with the spreading function /(Az) of one run, 

/ OO 

<t>{z)f{z- z)dz (16) 

-OO 


(where the particles which would penetrate the walls have 
to be “folded back” to the wall, i.e. (j)s{d) = P{z)dz). 
For constant run length, the spreading function /(Az) is 
obtained from the transfer functions as 


fiz) = r 

J z 


P{z'] 


dz' for z > 0, 


(17) 


and z < 0 follows from symmetry Note that if the tum¬ 
bling time Tt cannot be neglected compared to the run 
time Tp, the tumbling density has to be added propor¬ 


tional to T(/tp in Eq. (16). As an example, we consider 
here the case of thin three-dimensional channels {2d < L) 
and constant run length, which can be solved analytically. 
Here, the spreading function is found to be 


Note that the last two terms on the right-hand side of 
Eq. (15) determine the spatial dependence of the tumbling 
density near the wall. Figure (center) shows that this 
first-order calculation can qualitatively explain the numer- 


f{z) = 


2L 


In 


(18) 


^In the more general case, Eq. has to me modified to account 
for the run-length distribution 
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which yields the particle density (for z > 0) 


=4L 


In 

d 

4L2 


L 


d + 
2 + 


+ ln 


In 


d — z 
d — z 
d + z 


+ - d) + + d) 


(19) 


where the surface density ps of particles is obtained from 


the normalization condition. Equation (19) fits the simu¬ 


lations perfectly, without any adjustable parameters, see 

Fig-i 

Conclusions. — We have shown that run-and-tumble 
dynamics of self-propelled particles leads to highly struc¬ 
tured density distributions near impenetrable surfaces. 
Due to the absence of translational diffusion, accumula¬ 
tion materializes in the form of J-function peaks at the 
surface. Diffusion would broaden these peaks, similarly as 

Glose to confining walls, RTFs 


predicted for ABPs 26 


are thus clearly not equivalent to either diffusing particles 
or ABPs. The density distributions are predicted to de¬ 
pend sensitively on the spatial dimensionality and on the 
run-length distribution, where the typical length scale is 
set by the (average) run length. 

While the dynamics considered here is certainly over¬ 
simplified for real microswimmers like E. col% it captures 
the essential aspects of run-and-tumble motion, and simi¬ 
lar results can be expected for other types of Levy flights. 
In particular, the limit of large wall separations for the 
accumulation is very generic, and should thus apply to 
many systems of self-propelled particles and microswim- 
It will be interesting to see whether this 


mers 


21 22 


behavior extends to systems in which hydrodynamic in¬ 
teractions play a significant role. 

Another interesting issue is the behavior of RTFs at 
finite particle density [9 -11 . For high density, the char¬ 
acteristic features in confinement revealed by our study 
will almost certainly be washed out, because collisions will 
dominate over tumbling events. However, interesting be¬ 
havior can be expected in confinement, when the average 
distance between particles becomes comparable with the 
run length. 
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